Guiding-center dynamics of vortex dipoles in Bose-Einstein condensates 
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A quantized vortex dipole is the simplest vortex molecule, comprising two counter-circulating 
vortex lines in a superfluid. Although vortex dipoles are endemic in two-dimensional superfluids, the 
precise details of their dynamics have remained largely unexplored. We present here several striking 
observations of vortex dipoles in dilute-gas Bose-Einstein condensates, and develop a vortex-particle 
model that generates vortex line trajectories that are in good agreement with the experimental data. 
Interestingly, these diverse trajectories exhibit essentially identical quasi-periodic behavior, in which 
the vortex lines undergo stable epicyclic orbits. 



Introduction. — Vortices are persistent circulating flow 
patterns that occur in an extraordinary variety of sci- 
entific and mathematical contexts [1, 2]. In superfluids 
and superconductors, quantized vortices are topological 
excited collective states that play central roles in trans- 
port and dissipative properties. Dilute-gas Bose-Einstein 
condensates [3, 4] have supplied a particularly pristine 
setting for the study of static and dynamic vortex con- 
figurations, ranging from solitary vortex lines to large 
vortex lattices [5, 6]. Our primary focus in this Rapid 
Communication is on the behavior of the vortex dipole 
— a self-assembled, counter-circulating pair of quantized 
vortices — which lies at the heart of phenomena asso- 
ciated with quantum turbulence [7] and two-dimensional 
(2D) degenerate Bose gases [8]. Only recently have vortex 
dipoles been realized experimentally in dilute-gas Bose- 
Einstein condensates [9-11], and no unifying theoretical 
perspective on their dynamical behavior has yet emerged. 

Vortex dipoles in weakly-interacting superfiuid systems 
were first considered as one of several stationary con- 
figurations of vortex clusters occurring in harmonically 
confined, quasi-2D geometries [12, 13]. Subsequent theo- 
retical studies determined their energy and angular mo- 
mentum [14], and two independent modes were identified 
under small perturbations: a zero-frequency, precessional 
mode [15], and an oscillatory mode about the station- 
ary fixed points [16]. Additional stationary vortex clus- 
ter configurations, such as vortex tripoles [17] and other, 
more exotic arrangements [18] have also been predicted. 

Experiments involving vortex clusters are of more re- 
cent vintage [10, 11, 19]. Quantized vortex dipoles and 
the trajectories of their individual vortex lines were first 
observed experimentally by translating a repulsive obsta- 
cle through a Bose-Einstein condensate [10]; these obser- 
vations were accompanied by an integration of the three- 
dimensional (3D) Gross-Pitaevskn (GP) equation [3], 
which achieved good qualitative agreement with the ex- 
perimental trajectories of the individual vortex lines. The 



development of a real-time vortex imaging method led to 
the subsequent identification of asymmetric and station- 
ary symmetric vortex dipole configurations [11]. A recent 
theoretical analysis [20] of these results found good agree- 
ment with the static properties of the symmetric dipole, 
but met with less success in explaining the dynamics of 
the asymmetric dipoles when using the 3D GP model. 

We present here a detailed experimental study of vor- 
tex dipole dynamics in trapped Bose condensates, and 
we develop a theoretical framework that discerns the 
fundamental patterns of vortex motion by concentrat- 
ing directly on the vortices themselves. We find that, 
to lowest order, the vortices in a dipole behave as classi- 
cal particles, exhibiting quasi-periodic dynamics in which 
they execute stable, cyclic orbits about rigidly precess- 
ing guiding centers. Ultimately, we achieve an intuitive 
understanding of the underlying dynamics, in which the 
compelling agreement between theory and experimental 
results reveals that all such vortex dipoles exhibit essen- 
tially identical quasi-periodic behavior. 

Experiment. — Our apparatus and experimental pro- 
cedure are described elsewhere [11, 21]. We begin with 
an oblate ^^Rb condensate of ^ 6 x 10^ atoms in the 
\F — l,mF = —1) state, confined in a TOP magnetic 
trap [22] with effective frequencies {a;^,^^} ~ 27r • 
{35.8(2), 101.2(5)}Hz. The Thomas-Fermi radius of the 
condensate is R± ^17 /J,m. Vortex dipoles arise sponta- 
neously [23] but infrequently (~ 3% of the time) during 
the final stage of evaporative cooling. To observe the 
vortex dynamics we use a real-time imaging method [11] 
that repeatedly extracts, expands, and absorptively im- 
ages small (1% or 5%) fractions of the atomic sample. A 
final timc-of-fiight image of the remaining atoms is taken 
at the end of the experiment. 

Image sequences of the motion of the vortex lines are 
shown in Figs. 1 and 2. A surface fit to each image deter- 
mines the coordinates of each vortex and the condensate 
center. Vortex coordinates that cannot be fit are approx- 
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imatcd by direct inspection. We destructively measure 
the initial number of atoms in a typical condensate at 
the end of the evaporation cycle, and measure the final 
number of atoms in subsequent condensates that bear 
vortex dipoles. We estimate the number of atoms at in- 
termediate times based on the fraction extracted for each 
image, including an extra reduction (typically 1%) be- 
tween images to reflect atomic losses due to three-body 
recombinations [24]. 

The sudden reduction in the number of atoms follow- 
ing expansion is expected to excite a small-amplitude (a 
few percent) breathing mode, with only slight alteration 
of vortex trajectories [20]. Experimentally, we discern 
neither the anticipated excitation nor any disturbance of 
the motion of a solitary vortex line [11], further suggest- 
ing that any effect on a dipole is small. A full exploration 
of the imaging process and its effect on dipole dynamics 
awaits a theory of extraction and is therefore beyond the 
scope of the present work. 
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FIG. 1: (color) Symmetric vortex configuration, (a) Experi- 
mental images of the atomic density (optical depth), 1% ex- 
traction, (b) Polar plot of the experimental data (triangles) 
and the theoretical prediction from Eq. (2) (solid line). The 
colored circles correspond to the calculated vortex coordinates 
at the same times as the experimental data. Downward (up- 
ward) triangles indicate configurations in which a surface fit 
to the images in (a) succeeded (failed) to distinguish the vor- 
tices. Only the inner region of the condensate is shown. 

Theory. — We adopt an effective particle model in or- 
der to understand and theoretically follow the fundamen- 
tal characteristics of the observed vortex dynamics and 
interactions. This type of modeling has been shown to 
be successful in describing single vortex dynamics [5, 25] 
and, more recently, in capturing some multi-vortex fea- 
tures [16, 18, 26]. Particle modeling has also proven use- 
ful in identifying equilibrium states, oscillatory modes, 
and far-from-equilibrium collisional dynamics of matter- 
wave dark solitons [27-29]. 

The model for a vortex pair involves only two ingre- 
dients, and is readily generalized to n vortices [18]. The 



first ingredient is the gyroscopic precession of a solitary 
vortex line induced by the inhomogeneous atomic density 
profile of the condensate. In the absence of dissipation, 
this amounts to uniform circular motion at a fixed dis- 
tance r from the condensate center. For a nonrotating, 
axisymmetric, disk-shaped trap, the angular frequency 
51 (r) of the precession is given by [5, 11]: 
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where (w^) is the radial (axial) trap frequency, is 
the chemical potential, and R±^ is the Thomas-Fermi ra- 
dius. The second ingredient is the pairwise interaction 
between the vortex lines. For two vortex lines (labeled by 
subscripts j and k) in a homogeneous condensate, each 
point on one vortex line moves in the direction of the 
fluid flow at that point due to the other vortex line [30]. 
It is useful to express this interaction in terms of the 
frequency at which two co-rotating vortices orbit one an- 
other, $(^12) = ft/(mr^2)i where m is the atomic mass 
and ri2 is their separation. The linear speed of a dipole 
is u = ri2$(ri2). 

Each vortex in a "gas" of n quantized vortices moves 
in response to the sum of two velocities. 
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where Zk = Xk + iyk = r^e^^'' is the complex coordi- 
nate of the fcth vortex, Sk ~ ±1 is its topological charge, 
and rjk = \zk — Zj\. The constant parameter b modifies 
the interaction strength slightly from the homogeneous 
case; for our experimental conditions, b = 1.35 [18]. This 
lowest-order model does not take into account effects such 
as the bending of the vortex lines within the oblate con- 
densate or the dissipative interactions that occur between 
vortex lines and thermal atoms at finite temperature [31]. 

Analysis. — The many vortex dynamics described by 
Eq. (2) provide a number of experimentally accessible 
predictions for the observed symmetric and asymmetric 
vortex dipole configurations. First, an explicit prediction 
can be made about the stationary equilibrium positions 
of the constituent vortices. The velocities of the two vor- 
tices are zero for all times at a separation distance: 
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with flo = 51(0) and = ^{Ri_). Equation (3) yields an 
average equilibrium distance of 0.436i?j_, in remarkable 
agreement with the experimental value 0.432(17)i?_L [11]. 

The frequency of the epicyclic motion observed in 
Fig. 1 is calculated by linearizing the two-vortex equa- 
tions around dsep- We obtain two pairs of frequencies, 
a;i.2 = 0, and 
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FIG. 2: (color) Asymmetric vortex configurations, (a-c) Ex- 
perimental images of the atomic density overlaid with calcu- 
lated vortex positions from the model, Eq. (2) (colored rings). 
The images in (c) originally appeared in Ref. [11]. (d-f) Tra- 
jectories for the data in panels (a-c), using the notation of 
Fig. 1. (g-i) Trajectories for the data in panels (d-f) in the 
indicated rotating frames. In (a-b) the extracted fraction is 
1%, whereas in (c) it is 5% and the image intensity is scaled 
by a factor of two. 



The first pair illustrates the rotational invariance (around 
the center of the condensate) of the dipole [15], whereas 
the second pair corresponds to the anomalous modes 
of the vortex system [5], similar to those of dark soli- 
tons [27-29]. These anomalous modes are responsible 
here for the counter-rotating, oscillatory vortex motion 
about the equilibrium points [18], which can be under- 
stood intuitively as the result of competition between the 
precession and interaction velocities [16]. The measured 
oscillation frequency is 1.66(5) fio, where the overbar de- 
notes the average value during the experiment. This 



compares reasonably with the prediction of 1.52 fig from 
Eq. (4), given the relatively large amplitude of the ob- 
served vortex motion. The initially cyclic vortex motion 
observed in Ref. [10] is also of this persuasion, albeit with 
larger amplitude. 

A direct comparison between experimental trajectories 
and those predicted by the model is made possible by in- 
tegrating the equations of motion (as shown in Fig. lb). 
It is convenient to make this comparison by rescaling 
the post-expansion data by the ratio R± / R*^ , where R*^ 
is the measured condensate radius after expansion, and 
R± is the theoretical in-trap Thomas-Fermi radius cal- 
culated for a condensate of N{t) atoms. The theoretical 
points are calculated from Eq. (2), using values for the 
chemical potential that are also determined from N{t), 
and are updated after each measurement in the subse- 
quent time-evolution. We choose the initial positions of 
the modeled vortices to coincide with those of the experi- 
ment at t = 0. We also include, as a sole free parameter, a 
small (< 1 fim) offset (i.e., uncertainty) of the condensate 
center {Sx, 5y) with respect to all of the measured data 
points. The offset is attributed to a technical effect aris- 
ing from atomic micromotion in the time-averaged po- 
tential of the magnetic trap; similar shifts have been ob- 
served in condensates with single vortex lines [11]. With 
these parameters and initial conditions, we find good 
agreement between the experiment and the results of the 
numerical integration (Fig. lb) for symmetric displace- 
ments of the vortices from the equilibrium points. 

Next, we consider vortex trajectories arising from 
asymmetric displacements of the vortex lines. Three typ- 
ical image sequences are shown in Fig. 2a-f. In Fig. 2d, 
cyclic motion akin to that of Fig. 1 is superposed upon 
a slow precession. With increasing asymmetric displace- 
ments. Fig. 2e-f, one vortex moves in loops near the cen- 
ter of the condensate, while the other orbits closer to the 
periphery. The agreement between the simulation and 
the experiment ranges between fair (d) and quite good 
(e-f), with more pronounced discrepancies when the vor- 
tices are close together and difficult to resolve. 

Interestingly, all three data sets of Fig. 2 exhibit es- 
sentially identical epicyclic behavior, as is easily seen by 
transforming in each case to a suitable co-rotating refer- 
ence frame (Fig. 2g-i). We therefore arrive at a unified 
and intuitive general understanding of the vortex dipole 
behavior: the vortex lines oscillate in phase about sta- 
ble equilibrium points (or "guiding centers" ) that pre- 
cess rigidly about the center of the condensate. Both the 
precession frequency and the cyclic frequency (in the ro- 
tating frame) increase with displacement from the con- 
densate center, and the amplitude of the outer vortex 
orbit becomes smaller than that of the inner vortex; this 
last fact is due to the faster precession of the outer vortex 
chiefly influencing the co-rotating frame frequency. The 
stationary vortex dipole, as well as the oscillating dipole 
in Fig. 1, are special cases of this general motion in which 
guiding centers are stationary in the laboratory frame. 

Such guiding center motions can be mathemati- 
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cally identified in the context of the underlying model 
(Eq. (2)). Assuming rigidly rotating solutions zi{t) = 
riexp(iujt) and Z2{t) = — r2 cxp(ia;t), in terms of the 
fixed distances of the vortices from the condensate center 
ri and r2 (in units of R±), we find the following expres- 
sion for the precession frequency cj: 
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where bo = hb/{mRl), a = 1/(1 - rj), P = 1/(1 - rl) 
and 7 = (ri + r2)^^/2. A polynomial equation can also 
be derived for ri (and r2). A (nonlinear) analysis of such 
equilibria illustrates their stability, thus confirming the 
generic quasi-periodic nature of the vortex motion. 

Conclusions & Outlook. — We have shown that a sim- 
ple vortex-particle model reproduces the main features of 
vortex dipole behavior, as is experimentally observed in 
oblate Bose-Einstein condensates. The theoretical pre- 
diction of the separation distance between the station- 
ary fixed points and the frequency of small oscillations 
around these equilibria are in good agreement with the 



experimental data. We ultimately understand the generic 
vortex behavior and its quasi-periodic character as cyclic 
oscillations about a broader class of effective equilibria 
consisting of rigidly rotating guiding centers. Stationary 
or periodic motions arise as special cases within this set. 

Extensions of this simple vortex-particle approach may 
be applied to larger ensembles of vortices, including co- 
rotating vortices that have not crystallized into a lat- 
tice. By improving the ability of the apparatus to image 
the extracted atomic samples, we can contemplate the 
intriguing possibility of detecting additional vortex in- 
teractions through departures from the model, such as 
the effects of bent vortex lines, reconnections, and other 
phenomena associated with quantum turbulence. 
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